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Kafri, Mukamel, and Peliti Reply: The Comment 
of Hanke and Metzler jlj] questions the vaUdity of the 
analysis presented in Q to DNA chains of finite length as 
used in experiments. Their argument is that for the anal- 
ysis to be valid "each of the three segments going out of a 
vertex must be much longer than the persistence length 
(.p of this segment". By using the persistence lengths 
£p(L) ~ 40 A for a single strand and ip{H) ~ 500 A for 
a double helix (bound segment) they arrive at the con- 
clusion that in order to observe the asymptotic behavior 
found in 0] one needs chains which are far longer than 
those studied experimentally. 

This assertion is a consequence of a misinterpretation 
of the analysis given in |2] . In this analysis one considers 
a loop interacting with the rest of the chain and not just 
with the vicinal double helices. Thus, in (21 each of the 
two lines attached to the loop is in fact composed of an 
alternating sequence of bound segments and denaturated 
loops. It may be viewed as a stick-and-joint structure, 
whereby adjacent double helices are loosely attached to 
each other via an open loop. The "rest of the chain" as 
considered in reference is in fact a jointed rods struc- 
ture with a persistence length which is given by ip{H) 
when the distance between loops, is larger than £p{H) 
and is given by ^ when < ip{H). 

In the Comment the authors claim that the analysis 
of is valid only when the helical segment is flexible 
and thus much larger than the persistence length ip{H), 
where it could be described by a self-avoiding walk. This 
claim is incorrect. The analysis is valid as long as the 
chain contains a sufhcicnt number of loops to allow con- 
sidering the stick-and-joint structure (Fig. 1) as a self- 
avoiding walk. Thus, even short chains of length of the 
order of 5000 base pairs could allow for about 20 joints 
which is large enough to make the analysis of 2] valid. 
This is clearly evident from experimental melting curves 
0,01 of chains of such lengths. These curves exhibit over 
10 individual sub transitions corresponding to the exis- 
tence of a large number of loops (since some of the peaks 
correspond to more than one loop). Thus, the claim made 
in that the Fisher exponent c — 1.766 should be em- 
ployed when the bound segment is a rigid rod is not valid. 
Moreover, many experiments were carried out on chains 
of a length of the order of 10^ base pairs (Tjj where the 
asymptotic regime considered in is easily accessible. 

Our observation ^ that the effective loop entropy pa- 
rameter c is in fact larger than 2 affects other param- 
eters which have been applied in modeling DNA. One 
such example is the cooperativity parameter ao, which 
affects the melting curve rather drastically. Within the 
Poland-Scheraga (PS) approach the average distance be- 
tween loops near the transition is given by ^ ~ l/co for 
c > 1. This parameter has been estimated in various 
studies to be very small, of the order of lO^'^ to 10^'"' 


^^-i/(2-c)^ The large 
observed ^' yields a very small (Tq when fit with, for ex- 


0, 0, H, 0, ■ We point out that this small numerical 
value is obtained by fitting experimentally sharp melting 
curves to a theory which yields a continuous transition 
0,0 (namely, the PS model with 1 < c < 2), or some- 
times no transition at all (c = 0). Small ctq makes a 
continuous transition look sharp, yielding a good fit with 
experiment. In particular, ctq is usually estimated by con- 
sidering the maximal value of the temperature derivative 
of the order parameter 6 near the transition. For c < 2 
one has i * = - (9e/9r),„,, oc a' 
'ery sn 

ample, c = 1.5, 1.766 or c = 7]. However, since 

the effective c is larger than 2, the transition is expected 
to be first order where \E' is infinite, irrespective of ctq. 
The large but finite ^' values found in experiment could 
be attributed to finite size rounding effects of a first order 
transition rather than to a very small ctq. 

We believe that the vast amount of melting curves data 
existing in the literature should be reevaluated using the 
correct exponent c in order to obtain a realistic estimate 
for the cooperativity parameter. Certainly this should re- 
sult in a larger cooperativity parameter ctq which would 
in turn be reflected in an even smaller persistence length 
of the stick-and-joint structure of the molecule. It would 
also be very interesting to measure ctq directly from sin- 
gle molecule experiments using, for example, fluorescence 
correlation spectroscopy. 

We thank H. A. Scheraga and E. Yeramian for helpful 
discussions. 
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